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Abstract 

We study the CP violation effects from two types of neutrino mass matrices with (i) ( M v ) ee = 0, 
and (ii) ( M u ) ee = (M u ) efl = 0, which can be realized by the high dimensional lepton number 
violating operators Ll(D and £ c r Ir{D^<- h) 2 d> 2 , respectively. In (i), the neutrino mass 
spectrum is in the normal ordering with the lightest neutrino mass within the range 0.002 eV < 
m 0 < 0.007eV. Furthermore, for a given value of mo, there are two solutions for the two Majorana 
phases «2i and 0:31, whereas the Dirac phase 5 is arbitrary. For (ii), the parameters of mo, <5, 0:21 > 
and «3i can be completely determined. We calculate the CP violating asymmetries in neutrino- 
antineutrino oscillations for both mass textures of (i) and (ii), which are closely related to the CP 
violating Majorana phases. 
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I. INTRODUCTION 


Although it has been established that neutrinos are massive and mix each other in the 


recent several decades 


-7], their nature is still mysterious. It is known that neutrino mass 
terms could be of the Dirac type, in analogy to the charged fermions, i.e. quarks and charged 
leptons, or the Majorana type, possibly generated by the Weinberg operator L c L < h$ |8]. In 


the literature, there have been many models to realize the Weinberg operator at the tree 


0 - 


21| and loop 22428] levels. Note that the current neutrino oscillation experiments cannot 


determine the three CP violating phases, especially for the two Majorana phases, which is 
an important problem in neutrino physics. 

The Weinberg operator violates the lepton number symmetry by two units, but sometimes 
it is not the one that gives the dominant contribution to the Majorana neutrino masses 


rer dimensional LNV 


32], and dimension-9, 


or the lepton number violating (LNV) processes. Instead, other hi 
operators, for examples, the dimension-7, = £ R 'y fJj L L (D ti &)& 2 [29 
Og = I c r £r(D ^) 2 ^ 2 33-38], operators can lead to new Majorana neutrino mass structures 
different from those by the Weinberg operator if they are prominent. Specifically, due to 
the nontrivial dependence of the charged lepton masses, Oy generically generates a neutrino 


mass matrix with (M u ) ee = 0 [32|] in the 
texture with (M u ) ee = (M„) e/t = 0 [28, 29 


iavor basis, while Og naturally gives rise to the 
33L 1351 . 1371 [38]. By fitting the present neutrino 


oscillation data, both textures predict that the neutrino mass matrix should be of the normal 
ordering, and already give stringent constraints to the unknown parameters in the neutrino 

I, we show that they could naturally lead to 


mass matrix. In particular, in Refs. 32 


nontrivial values for the three CP violating phases. We regard that the higher dimensional 
operators would provide us with a new way to generate the new neutrino structures, besides 


the ordinary approach by imposing flavor symmetries 


39 


49]. 


In the present paper, we investigate a relevant question: to what extent can the conditions 
(M u ) ee = 0 and (. M v ) ee = (M u ) efl = 0 restrict the neutrino mass matrix structure, especially 
the leptonic CP violating phases, based on measured quantities from oscillation experiments? 
In our treatment, we also take into account the experimental uncertainties in the data 
in order to see their effects on the results. Note that there were already many studies 
about these two specific neutrino mass matrices in the literature, see e.g., Refs. 50-67], in 
which more texture-zero neutrino mass matrices were examined. Here, our focus is their 
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implications on the leptonic CP violating phases. 


With the predicted neutrino mass parameters, the next question is how to test the above 
two texture-zero structures by measuring all the relevant parameters in the neutrino mass 
matrices, especially the non-trivial Majorana phases. Previous studies showed that neutrino- 


antineutrino oscillations gave us prospective approaches to probe the Majorana phases 


68 - 


77], which is impossible for the conventional (anti)neutrino-(anti)neutrino oscillation exper¬ 


iments. We find that once the possible regions of these phases are depicted for the present 
two textures, the associated CP violating asymmetries of the neutrino-antineutrino oscilla¬ 
tions can be predicted. As will be shown later, by appropriately choosing the (anti)neutrino 
beam energy and baseline length, some of the asymmetries can be of 0(1). 


This paper is organized as follows. In Sec. II we study the implications of the texture- 
zero conditions (M„) ee = 0 and (M„) ee = (M v ) efl = 0 to the unknown neutrino mass 
parameters, including the lightest neutrino mass and the three CP violating phases, based 
on the existing data. With the preferred values of these parameters, we predict the CP 
violating asymmetries in the neutrino-antineutrino oscillations for both textures in Sec. III. 
In Sec. IV, we give a short summary. 


II. TEXTURE-ZERO NEUTRINO MASS MATRIX 


As the Majorana neutrino mass matrix is symmetric, there are six independent com¬ 
plex elements (M„) ee , (M,,)^, and (. M V ) TT . A well-defined M v can 

be connected with the observed quantities from neutrino oscillations. Up to the held redef¬ 
inition, all of the above matrix elements depend on the nine neutrino parameters, including 
3 masses, 3 mixing angles, 1 Dirac CP phase and 2 Majorana CP phases. In the flavor basis, 
where the charged lepton mass matrix is diagonal, the neutrino mass matrix defined in the 
Lagrangian L = — \{^u) c Mw vjjj + H.c. can be decomposed as follows, 


M u = V*diag(mi, m 2 , m 3 )U t , 


( 1 ) 
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where 777 . 1 , 2,3 are three neutrino masses. V is the charged current lcptonic mixing matrix 


78 


79], conventionally expressed in the standard parametrization as 80.1811 



C12C13 S12C13 s 13 e~ tS 


h ^ 

O 

O 

V = 

~ S12C23 — Ci 2 S23 s 13 e * <5 D2C23 — -Si2S23Si3e i<5 S23 c 13 


0 e ia21 ' 2 0 


1 S12S23 — Cl2C23 s 13 e4<5 — C12S23 — S 12 C23S 13 e tS C 23 C 13 , 


0 0 e ia 3l/ 2 


, ( 2 ) 


where s l3 = sin 9,j , c„- = cos 6V, S is the Dirac phase, and « 2 i, 3 i represent two Ma- 


jorana phases within the range [0, 2tt\. The values of $ 12 , 0 23 , $ 13 , Am^ and |A?77 


have already been obtained from the neutrino oscillation experiments 
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81], so that the 


rest four unknown neutrino parameters are the three CP phases, 5, a 2 i and 032 , and 
the lightest neutrino mass, m 0 . Note that only the absolute value Am | 2 has been ac¬ 
quired, which leaves us two possible orderings: the normal ordering for Am | 2 > 0 with 
( 7771 , 7772 , 7773 ) = (mo, + A? 772 i, \/ m o + A/ 77 ^), and the inverted one for Am^ 2 < 0 with 
(7771,7772,7773) = (a/t 77 o - A77733, V<^^L + ^ fi , rno ). 

M u can have some special approximate texture-zero forms when it is generated by some 
high dimensional LNV operators. For example, if O 7 = £ c r 'J^Ll(D i 1 Q)Q 2 gives the leading 
contribution to neutrino masses, then ( M v )w should be approximately proportional to the 
sum of charged lepton masses, mi + m^, with £ and i' = e, /i, r. Consequently, (M u ) ee 
should be much smaller than other elements. Similarly, if Og = £ r £r(D^) 2 Q 2 dominates 
over other LNV operators, (M u )ui will be proportional to mimi'. It turns out that not 
only (M u ) ee but also (M u ) efl are expected be greatly suppressed due to the hierarchy in the 
charged lepton masses. In other words, the neutrino mass matrices obtained from these 
LNV effective operators are characterized by the special zero textures (M„) ee = 0 and 


(M u ) ee = (M u ) eii = 
two subsections. 


a 


the implications of which will be discussed in detail in the following 


1 Besides the relative smallness of the element ( M„) ee already argued in the main text for the two high¬ 
dimensional effective operators of O 7 and Og, its absolute value is further constrained by the neutrinoless 


double beta (Ov/3/3) decay processes 


82 


881 ]. Note that these two effective operators give the dominant 


contributions to the (V/3/3 decay at tree level, while the Majorana mass terms arising from CV( 9 ) begins at 
one-(two-)loop level. Due to the absence of the loop suppression, these two operators are more sensitive 
to the Oi'd/J decay processes, which constrain the cutoff scales and Wilson coefficients of the effective 


operators greatly and lead to the conclusion that (M„) ee < 10 13 eV. For further details, please refer to 
and 38]. 


operaiors 
Refs. 3 
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A. {M v ) ee = 0 


By expanding the right-hand side of Eq. ([I]) with the standard parametrization of V in 
Eq. (|2]1 . the condition (M u ) ee = 0 can be transformed into the following relation, 


(M u ) ee = c\ 2 c{ 3 mi + s‘i 2 ci 3 m 2 e 10121 + s{ 3 m 3 e lA = 0 , 


( 3 ) 


where the phase A = a 3 1 — 25 is defined, which will be used to replace ^31 as an independent 
Majorana phase hereafter. Note that this equation excludes the inverted ordering at more 
than 2 a significance by current oscillation experiment results [81], so that we only need to 


consider the normal-ordering neutrino mass matrix from now on. Fig. Q] shows the allowed 
parameter space region satisfying Eq. (j3J), in which the solid curves represent the parameters 
when the experimental observables are at their central values in PDG 81], while the shadow 
areas correspond to the la standard deviations, ft is interesting to note that mo and a 2 i are 
already limited within small parameter regions, with 0 . 87 T < a 2 i < 1.27r and 0.0015 eV < 
m 0 0.008 eV, respectively. In particular, the extremal values of the lightest neutrino 
mass (rno)min(max) are related to two CP conserving solutions to Eq. ([3]) with a 2 \ = 7 r and 
A = 0 (tt), 


s 2 i 2 c 2 i 3 m 2 - c\ 2 c\ 3 m\ - s\ 3 m 3 


0 and c 2 l2 c\ 3 rrii — s\ 2 c\ 3 m 2 — s\ 3 m 3 


0 . 


( 4 ) 


For each value of m 0 within the regions [(mo)min, (wo)max], there exist two solutions for a 2 \ 
and A, differentiated by the positive or negative sin am, which are shown in Fig. |T] as red 
or blue curves/shadows. Another interesting observation is that the obtained 021 is limited 
around 7 r, which can be understood directly from Eq. ([3]). Since s 33 is very small, the third 
term in Eq. ([3)) can be neglected, and the first two terms must balance each other to achieve 
the constraint of the vanishing (M u ) ee , which only requires a 2 i ~ n in order to reverse the 
sign of the second term. Moreover, a 2 i is precisely predicted to be 1.17T or 0 . 97 T when mo is 
located within 0.004 eV < mo 0.005 e\u, no matter how experimental errors vary. Finally, 


we remark that ( M u ) ee = 0 does not provide any constraint on 5, which is only contained in 
A. If one focus on the real M u , then 5 can be taken as 0 or 7 r, for the cases mo — (^cOmm and 
(mo)max- Therefore, there are 4 independent real neutrino mass matrices for (M„) ee = 0. 


2 Similar results are also given in Ref. 
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FIG. 1. Contours for (M u ) ee = 0 in (a) mo-a^i and (b) mo-A planes, respectively. Solid lines with 
experimental central values and la standard deviation, where the colors of red and blue indicate 
that the values of sin « 2 i are positive and negative, respectively. 


B. (M„)ee = = 0 


In this subsection, we will concentrate on the case with (M„) ee = (M v ) efl = 0. Note that 
such constraints correspond to two complex equations, which enable us to uniquely solve 
for the remaining four parameters (mo, S, a 21 , A) in the neutrino mass matrix undetermined 
from the current oscillation experiments. Now we sketch the procedure of deriving these 


quantities in terms of the measured observables 


28]. The first step is to write down the two 


conditions in the parametrization independent form: 

{M u )ee = m-lVel + = 0 , 

(M u ) e , = mi v e to + m 2 v; 2 v ; 2 + m 3 v; 3 v ; 3 = o, 

with which we can obtain the following useful formulas 

1 - IH 2 


Amli 


1 - |X| 2 ’ 

m 2 |y| 
m 3 \X\ ’ 

im(v e 3 v r 3 v; 2 v; 2 ) = im(v el v rl v; 3 v; 3 ) = im{v e 2 v T 2 v: x v T \) = o, 


(5) 


( 6 ) 

(7) 

( 8 ) 


with 


X = 


K 2 K2KlK*l 


Y- 1 = 


[VetflVeiV^ - V e3 V^\ 


v el v rl v; 3 v r3 


\v e 3 \ 2 \v e 3 v tt 2 -v e 2 v^\ 


(9) 
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Since neither |X| nor |Y| depends on the two Majorana phases and mo, we can determine 
the Dirac phase 5 from Eq. ()6|) . By substituting the obtained Dirac phase into Eq. o , we 


can solve for m 0 . Finally, two Majorana phases can be fixed with Eq. (jSl) . In the standard 


parametrization, the solution is expressed by 


28] 


cos S 
A 
a 2i 


m 0 


_ 'hi _ 

2(l + t? 2 )+r(l-f? 2 )U 23 113 
arg(-si 3 + t 12 t 23 e~ lS ) , 


/ gi 3 ~ t 12 t 23 e iS \ 
' S13 + ti 2 t 23 e~ lS y 


) 


i) 



(10a) 

(10b) 


(10c) 



(2 + r) tf 3 (tj 2 t 23 — 2 s 13 t 12 t 23 C 5 + sf 3 ) 

2 r 1 - tf 3 (l + tf 2 tl 3 - 2si 3 ti 2 t 23 c 5 ) ’ 


(lOd) 


with r = Am| 1 /(Am| 2 + A/n 21 /2). Note that for each value of mo, we can obtain two 
solutions of the CP violating phases (5, a 2 1 , and A), which can be connected with each 
other by the replacements of <5 — y 2n — 6, a 2 i — > 2tt — a 2 i, and A — > 2% — A. 


Fig. [2] shows the allowed parameter space according to Eq. (fTOl) when we take the fitting 
values of the five parameters (0 12 , 9 23 , 9 1 3 , Am| l7 Am 32 ) within the la deviation as the input 
parameters, in which the red (blue) areas label the regions with sina^i > 0(< 0). If we 
take the central experimental values of the measured quantities 81], two solutions can 
be obtained with m 0 = 5.07 x 10~ 3 eV, 5 = 0.597T (1.4l7r), cx 2 i = 0.897T (l.ll7r), and 
A = 1.347T (0.667r). The corresponding leptonic Jarlskog invariant, J = Ci 2 c^ 3 c 23 Si 2 Si 3 s 23 S 5 , 
is equal to 0.033(—0.033), which characterizes CP violation in the lepton sector. It is worth 
noting that the value of 5 = lAln for one of the solution is close to the central value of 5 
from the global fitting result 8J-] - Moreover, in the present texture-zero case, when Dirac 
phase 6 is taken to be the CP conserving values, such as 5 = 0 or 7r, the two Majorana 
ones can only be CP conserving values too, i.e., a 2 \ and A should be 0 or 7r. However, with 
the results shown in Fig. [21 it is interesting that the CP conserving cases are excluded at 
the la level. Finally, if the Dirac phase is taken to be of the maximal CP violating value 
with S = 7t/2 (— 7t/2), the Majorana phases are predicted to be ct 2 i = 0.887T (1.127t) and 
A = 1.427r (0.587t) with the experimental central values for the mixing angles. 
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FIG. 2. Correlations among the parameters: <5, A, 021 , and mo, where the red (blue) color 
represents that sin 0:21 > 0(< 0). 


III. NEUTRINO-ANTINEUTRINO OSCILLATIONS 


With the non-trivial Majorana CP violating phases, especially for the case with M ee = 
M efl = 0, the immediate important question is how to measure them. Traditional (anti)neutrino- 
(anti)neutrino oscillation experiments can be used to measure the Dirac phase, but they are 
insensitive to the Majorana ones since the involved processes are lepton number conserving 
so that the Majorana phases are cancelled out in the corresponding formulas. As a result, in 


order to measure the Majorana phases, one o: 
processes are LNV. As pointed out in Refs. 


le necessary conditions is that the involved 
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.71], the neutrino-antineutrino oscillations 


provide us with a promising way to detect them. Unfortunately, the neutrino-antineutrino 
channels suffer from an additional helicity suppression factor ( m v /E ) 2 in the oscillation 
probabilities, compared with the corresponding usual (anti)neutrino-(anti)neutrino oscilla¬ 
tion channels. Therefore, it is challenging to carry out such experiments. The use of the 
















low-energy Mossbauer electron antineutrinos 89] with 18.6 keV, which are emitted from 
the bound-state beta decay of 3 H to 3 He, can improve the situation greatly by enhancing 
the signal by a factor of 0 ( 10 4 ) as compared with the conventional reactor antineutri¬ 
nos [71]. However, even in this case, it is still practically impossible to observe these 
neutrino-antineutrino oscillations, as will be shown below. 

The general formulas for the neutrino-antineutrino oscillation probabilities P(u a —> Dp) 
and P(uf 3 —> D a ) in the three-flavor framework are 71] 


\K 


P(v a vp) = \M a p \ 2 - 4 ^ rn l rri J R,e(V ai Vg l V* J Vp j ) sin 2 ( 


i<j 


/Am^L 

4 E 


+2 ^ r n l m : jlm(V ai V (il V* rj ) sin ( 


/A m^L 


i<j 
12 


2 E 


( 11 ) 


P(v a -)• up) = S- [\M a p \ 2 - 4 m^ReiVarVpiV:^) si 

i<j 


sm 


/ ArrijjL 
V 4 E 


-2 ^ riknij Im( V m Vp t V*j ) sin 


i<j 


/Am 2 i L\^ 

2 E 


( 12 ) 


where K and K are the kinetic factors with \K\ = \K\ and L is the neutrino traveling 
length. Now it is interesting to estimate the neutrino-antineutrino oscillation probabilities 


for different channels to see if they have the potential to be observed unc 
experimental status, especially the Mossbauer neutrinos advertised in Ref. 


er the present 


71]. By assuming 


the kinematic factor K ~ 0(1), electron antineutrino energy E ~ 18.6 keV, and oscillation 
baseline length L ~ 300 m, we can obtain the largest u e — D e oscillation probability to be 
P(u e -> D e ) rs j 0( 10 13 ) for m 0 = 0.0065 eV. The largest probabilities for other oscillation 
channels, such as P(D e —> u^), would be of the similar order. In the view of these simple 
exercises, it seems impossible to observe these oscillations practically in the foreseeable 
experiments. 

It is obvious that P(u a —* Dp) and its CP conjugate process P(u a —» Dp) can have different 
values when V is complex, which is the origin of CP violation in the lepton sector. Therefore, 
we can define the CP asymmetry parameter A a p by 

_ P(u a ->• Dp) - P(D a ->■ up) 


A-ap — 


P(u a —> Up) + P(u a —> Up) 

2 E i<j mimjlmiVaiVpiV*^) sin ((fnAm 2 ji )/( 2 Am 2 21 )) 

\ M a p\ 2 ~ 4 J2i<j mmjReiVaiVpiV^V^) sin 2 ((/7rAm 2 J/(4Am 2 1 )) ’ 


(13) 
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where / = (L/F)(Am 21 /7r). In the following, we shall use the obtained CP violating 
phases from the previous two texture-zero neutrino mass matrices to predict the oscillation 
probabilities and the associated CP violating asymmetries in some neutrino-antineutrino 
oscillation channels of great phenomenological interest, and then see how these measurements 
can help us to probe or constrain the whole picture of neutrino masses. 


A. (M,,) ee = 0 


According to Eq. (TT3T) . the asymmetry A ee can be expressed by 
4U = 2 


m\rri 2 C 2 V2 c 4 ]:i Sy 2 sin 2 sin a 2 1 + rn l m :i Cy 2 Cy S s 2 vs sin 2 


LT"( 


Amir. 

1 H-— 

Am 2 21 


sin A 


9 9 9 . 9 f 7T A 777-00 . , * \ 

-m 2 m 3 c 13 s l2 s l3 sm 2Am2 ~ sm («2i - A) 




+4mim 2 Ci 2 Ci 3 s 32 sin 2 cosc^i + 4mim3C 2 2 c 2 3 s 2 3 sin 2 
+4m 2 m 3 c 2 3S 2 2 s 2 3 sin 2 ^^™ 2 2 cos(a 2 i - A) 




/tt A Amj 2 
4 V Am 21 


cos A 


(14) 


When imposing the condition ( M u ) ee = 0, both a 2 \ and A can be expressed as the functions 
of mo- Therefore, by Exing the factor / to some definite value, the CP violating asymmetries 
of various channels can also have definite values for every tuq. As an illustration, Fig. [3b 
gives the correlation between A ee and mo when / = 0.55, where we only take the central 
values of the measured quantities in our calculation. By comparing Fig. [3b and Fig. |T}i. we 
see that the detection of A ee directly implies the existence of a nonzero sino; 2 i, and their 
signs are positively correlated. Our direct calculation confirms this observation. We also plot 
the variation of A ee against the factor / in Fig. [3b by taking mo = 0.004 eV, which shows 
that if we can fine tune the beam energy E or the baseline length L to make an appropriate 
value of /, a large CP violating asymmetry in the 14-14 channel |A ee | ~ 1 can be obtained. 
Furthermore, note that the dependence of A ee on the Dirac phase 5 is only through the 
combination of A = 031 — 25, so that even if 5 vanishes, there can still be quite sizable CP- 
violating effects in the u e -P e oscillation experiment due to the compensation from <341 in A, 
which starkly shows the significance of the Majorana phases in generating the CP-violating 
effects. Finally, we make the estimation of the asymmetries in other neutrino-antineutrino 
oscillation channels, by simply taking / = 0.55, mo = 0.004 eV, and 6 = 0 or 7t/2, together 
with other measured observables at their central values, with the results shown in Table H[ 
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TABLE I. Asymmetries for neutrino-antineutrino oscillations with / = 0.55 and uiq = 0.004eV. 


u -n-ee -fie/A -^-er -^ L fir 

0 l.O(-l.O) 0.38(-0.38) -0.71(0.71) 0.19(-0.19) -0.21(0.21) 0.17(-0.17) 
vr/2 l.O(-l.O) 0.61(0.20) -0.20(-0.47) -0.26(0.25) 0.18(-0.19) -0.15(0.16) 




;« 0 (eV) 


(a) 


(b) 


FIG. 3. A ee as functions of (a) mo and (b) /, where the red (blue) corresponds to sina 2 i > 0 (< 0). 

B. [M v ) ee = (M u ) ei j = 0 

It was shown previously that all of the mass and mixing parameters can be fixed when 
(Mj,) ee = = 0. It follows that all of the CP violating asymmetries in the neutrino- 

antineutrino oscillations can also be determined. We remark that for this case A er always 
vanishes since each term in the summation of Eq. (TT3l) is zero as the consequence of Eq. (El) • 
Using the central values of neutrino mixing parameters from neutrino oscillations with / = 
3.5, we can predict A ee = 0.92, A efl = —0.09, A^ = —0.15, A /IT = 0.1, and A TT = —0.09. 
Therefore, the ry-Te oscillation is the most prospective channel to probe this neutrino mass 
texture. 


IV. SUMMARY 

We have studied the CP violating asymmetries and related LNV processes such as the 
neutrino-antineutrino oscillations under two types of the neutrino mass textures, ( M u ) ee = 
0 and ( M u ) ee = (M u ) efl = 0, realized by the high dimensional lepton number violating 
operators. For (M u ) ee = 0, there are two solutions of a 2 i and A for each value of m 0 , with 
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FIG. 4. Leading-order Feynman diagrams for rare LNV B meson decays induced by (a) O 9 and 
(b) O r . 


0.002 eV < m 0 < 0.007eV and an arbitrary value of 5. For (. M v ) ee = = 0, two 

solutions for free parameters (mo, S, a 2 i, A) can be obtained, in which one of them with 
6 = 1.417T is close to the global fitting result. The effect of the nonzero values of the two 
Majorana phases can be reflected by the related CP violating asymmetry parameters A a p 
in neutrino-antineutrino oscillations. I 11 the texture (M„) ee = 0, we find that a non-zero A ee 
can be obtained even if the Dirac phase <5 is switched off, and its sign is positively correlated 
to that of sina 2 i- For (. M v ) ee = = 0, a large values of A ee is predicted, while A er is 

always zero. 


It is interesting to consider other probes to the Majorana character of the neutrino masses, 
such as rare LNV meson decays. It is well-known that ordinary channels with Majorana 
neutrino mass insertions are too small to be observed in the near future. However, it is 
remarkable that the effective operators, such as 0 7 and O 9 , would give new leading-order 
contributions. For concreteness, let us consider the process B + — y n _ /i + p + . If Majorana 
neutrino masses are induced by Og, the dominant channel to this process is given by the 
Feynman diagram in Fig. [4^, as this tree-level diagram does not involve the tiny Majorana 
neutrino masses which arise at two-loop level via Og. However, with the model parameters 
fixed by the observed neutrino masses as in Ref. 38], a simple estimation shows that the 
typical branching ratio for this process is to be of 0( 10~ 25 ). Other LNV rare meson decays, 
like K + — > 7 r _ /i + /i + , would have even smaller branching ratios. Similar results can also be 
obtained for 0 7 from Fig. [4}d. As a result, it seems also impossible to measure such LNV 
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meson decays practically. 
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